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We investigate a two-way cascaded cavity QED system consisting of microtoroidal resonators 
coupled through an optical fiber. Each microtoroidal cavity supports two counter-propagating 
whispering-gallery modes coupled to single atoms through their evanescent fields. We focus on a 
pair of atom-microtoroid systems and compute the spectrum of spontaneous emission into the fiber 
with one atom initially excited. Explicit results are presented for strong-coupling and bad-cavity 
regimes, where the latter allows the effective atom-atom interaction to be controlled through the 
atom-cavity coupling and detuning: the atoms exhibit either collective spontaneous emission with 
no dipole-dipole interaction or a (coherent) dipole-dipole interaction and independent (single-atom) 
emission. This capacity for switching the character of the interaction is a feature of bi-directional 
coupling and connects our two-way cascaded system to work on one-dimensional waveguides. Build¬ 
ing upon our bad-cavity results, we generalize to many atom-microtoroid systems coupled through 
an optical fiber. 

PACS numbers: 42.50.Pq, 42.50.Nn, 03.67.-a 


I. INTRODUCTION 

With a view to realizing distributed quantum networks 
in which atomic or solid state qubits are coupled through 
traveling light fields over substantial distances, a number 
of new architectures for cavity quantum electrodynamics 
(CQED) are proposed These include monolithic 

microtoroids and microdisks, whose whispering gallery 
modes may couple to free atoms via their evanescent 
fields m , or quantum dots embedded in the resonators 
themselves |10l - [l^ . Highly efficient input and output is 
achieved by overlapping the evanescent fields of tapered 
optical fibers with the resonator modes. Fibers are also 
the natural choice for linking resonators, i.e., as a channel 
of communication between different nodes of a network 
of such microresonator CQED systems. 

Quite apart from applications to quantum information 
science, such systems offer an intriguing new environment 
in which to study basic CQED phenomena, e.g., modified 
atomic spontaneous emission, where distant atoms inter¬ 
act through the resonator and fiber fields. New features 
compared with conventional single-mode CQED include 
counterpropagating resonator and fiber modes, and the 
open-system nature of the fiber-mediated coupling, where 
the fiber acts as a common, broadband reservoir to the 
modes of the microtoroids. A natural modeling of this 
coupling is provided within the framework of cascaded 
open quantum systems [3, [3- Starting from such a 
model, we focus for the most part in this paper on a pair 
of two-way cascaded atom-microtoroid systems [l^— the 
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many-site extension is briefly considered—and compute 
the spectrum of spontaneous emission when one atom 
is initially excited. After deriving general expressions 
for all system outputs, including side emission from the 
atoms and scattering in the microtoroids, we concentrate 
on the readily accessible spontaneous emission through 
the counter-propagating hber modes. 

Our model, most generally, describes a rich dynamic; 
the initial excitation is shared between two atoms and 
four microtoroid modes. We discuss two limiting cases in 
detail: the strong-coupling regime, with the atom-cavity 
coupling constants much larger than the system decay 
rates, and the bad-cavity (so-called ‘Purcell’) regime, 
where the cavity decay rates are sufficiently large that 
an adiabatic elimination of the microtoroid fields may be 
performed. In the strong coupling regime, depending on 
the relative coupling phase, we show that either a vacuum 
Rabi oscillation develops between delocalized—shared 
between the two atom-microtoroid systems—atomic and 
photonic excitations, or localized atom-photon dressed 
states couple through the cavity decay into the fiber. 

In the bad-cavity regime, after adiabatic elimination, 
again depending on relative coupling phase (alternatively, 
the atom-cavity detuning), the dynamics simplify to yield 
strongly contrasting forms of the effective—mediated by 
the microtoroid and fiber fields—atom-atom interaction: 
the atoms either exhibit collective spontaneous emission 
through the fiber with no dipole-dipole interaction, or 
a fiber-mediated dipole-dipole interaction with cavity- 
enhanced single-atom spontaneous emission. We recover 
behavior seen in the interaction of distant atoms through 
a one-dimensional waveguide [l7l - [^ . 

We present our cascaded open quantum system model 
for a pair of atom-microtoroid systems in Sec. im and 
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derive general expressions for the spontaneous emission 
spectra at its various outputs in Sec. lIIIl We then discuss 
the contrasting behaviors in the strong-coupling regime, 
Sec. lIVi and the bad-cavity regime, Sec. El In Sec. El 
we derive a master equation for the atoms alone, which 
provides a compact and transparent description of the 
atomic dynamics, and connects our two-way cascaded 
system, in the bad-cavity regime, to a one-dimensional 
waveguide. The master equation is generalized to many 
two-way cascaded atom-microtoroid systems in Sec. I VII 
where similar—but not the same—contrasting behaviors 
as a function of coupling phase are discussed. Spectra 
for iV = 3 two-way cascaded atom-microtoroid systems 
illutrate these results. Sec. IVIII provides a summary and 
conclusions. 


II. THEORETICAL MODEL 


rate A pair of identical two-level atoms, transition 

frequency uja and spontaneous emission rate 7 a, couple 
to the WGM’s through evanescent fields with coupling 
constants gi, i = 1 , 2 . 


Atom Atom 



FIG. 1: Schematic of the cascaded atom-microtoroid system. 
Counter-propagating microtoroid modes couple through their 
evanescent fields to two two-state atoms and the input-output 
fields of an optical fiber. 


We begin by investigating a system of two microtoroids 
coupled through an optical fiber, as shown in Fig. [l] 
The microtoroids are distinguished by a label i = 1,2. 
Toroid i supports a pair of counter-propagating whisper¬ 
ing gallery modes (WGM’s) of frequency wc,, with each 
coupled to the fiber through evanescent fields at loss rate 
2Kex ■ Light is lost out the sides of the nhcrotoroids—e.g., 
due to scattering from imperfections—at an intrinsic loss 


The system is modeled as a two-way cascaded system. 
We introduce creation and annihilation operators Uj, bi 
(clockwise, counterclockwise) and aj, bj for the WGM’s, 
and raising and lowering operators af for the two-level 
atoms. Then with Lindblad operator D[0\- = 20 ■ — 

0^0 ■ — ■ 0^0, the master equatiom in a frame rotating 
at frequency uja, is written as |14. Ha. 


p — — i [Hi + H 2 + Hi2, p] + — {D[Ja]p + D[Ji)]p) 

+ f^ltHD[ai]p + D[bi\p) + K^\D[a2]p + D[b2]p) + ^ {D[ai]p + D[a2]p) , ( 1 ) 

where p is the reduced density operator of the microtoroid-plus-atoms system, and (setting h = 1) the Hamiltonian 
is defined by 


Hi = Ai(a|ai -I- fojfei) -I- {hia\bi -I- h*ib\ai) -|- {gla\(7^ + gicr^ai) + {gib\(j^ + g^a^bi), 
H2 = A2(a2a2 -|- 62^2) + (^202^2 + ^2^2^2) + (52®2'^2 92<^2 (32^2'^2 + 52*^^^2): 

with Ai = uJCi — UJA, and 

Hi 2 = i\J Kex kIx (e“*'^“aja2 - e*‘^“a|ai) -I- i\j k^x ^ix {e~‘"’^'‘b\bi - e*‘^'’6|&2)- 


( 2 a) 

( 2 b) 


( 2 c) 


Hi describes the first atom-microtoroid system, H 2 the 
second, and H 12 the coupling between them via the fiber. 
Within each microtoroid, the counter-propagating modes 
can couple directly with strength hi, and (jja,b are phases 
that account for the propagation distance between the 


microtoroids. The operators 


Ja = ^ 

\J 2k^-} aie 

\J 2Kex 0-2 ; 

(3a) 

Jb = ^ 

j2K^^b2 + e-^^^\ 

/ 2 /cix 61 , 

(3b) 


describe a two-way cascaded coupling between the modes 
of the two microtoroids mediated by the fiber, which is 
regarded to provide a common (broadband) reservoir for 
these modes. Note that with the propagation directions 
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coupled—through scattering between counterpropagat- 
ing modes (/ii,2) and from the atoms (31,2)—the model 
as written is an approximation, in so far as it neglects re¬ 
tardation associated with the time of flight between the 
microtoroids. The approximation is valid provided the 
flight time is much smaller than the timescale on which 
the system state undergoes significant change , which 
in practice puts a limit on the length of fiber connecting 
the microtoroids. 

To summarize, the model incorporates three different 
types of damping: loss from the WGM’s into the fiber 
at rate 2Kex, intrinsic scattering loss from the WGM’s 
at rate 2Kin, and atomic spontaneous emission to free 
space at rate 7 a. In addition, the connterpropagating 
modes of each microtoroid are directly coupled (hi 2) and 
each atom couples at equal strength to both modes of its 
respective microtoroid ( 91 , 2 )■ 


III. SPECTRUM OF SPONTANEOUS 
EMISSION 

A. Quantum trajectory unraveling 

In this section we wish to calculate the spectrum of 
emitted photons starting from an initial state with one 
atom excited and the other in the ground state. We use 
the method introduced in Chap. 13 of [2^, which writes 
the master equation as p = Cp and separates C into 
two parts: the first part acts only on the one-quantum 
subspace, and the second is the generator of transitions 
out of this subspace to the ground state. Effectively we 
implement a quantum trajectory unraveling and compute 
the spectrum from the response of a frequency selective 
detector placed in the output field (Sec. 19.3.2 of 0). 

With C decomposed in this way, the master equation 
is written as 

p={C + V)p, (4) 

with superoperator C acting only in the one-quantum 
subspace: 

C . = - i[H, +H 2 + ■ ] - llJlJa + 4 Jb, • ] + 

- + b\bi,-]+ - 4n^[4«2 +&J& 2 ,-] + 

+^2(^2^] + ^ (5) 

with [•, •]_!_ the anti-commutator, and superoperator V 

generating transitions to the ground state: 

V-= Ja-Ji + Jb-Jl 

-\- 2 K; 4 (ai ■ a\+bi ■ 6 j) -|- 2K-^\a2 ■ 4 -I- 62 ■ 4 ) 
+ -yA(cr^ ■ cr++a2 ■ a+). (6) 


The initial density operator with the first atom excited 
is written p(0) = PA(0)pa&(0)i with atomic state 

PA(0) = (|e)(e|)A,(|g)(g|)A., (7) 

and Pab(O) the vacuum state of the field. 

Now for times t > 0, a photon has either been emitted, 
probability T’spon(i), in which case the system is in the 
ground state, 

|G) = |g)Ajg)A.|0),j0)bj0)a,|0),„ (8) 

or no photon has been emitted and the system is in a 
pure one-quantum state: 

l'0(i)) = + ai{t)a\ + ^i{t)b\ 

+ ^2it)o-2 + a2t)al + P2t)bl]\G), (9) 

where ^i{t) is the probability amplitude for the first atom 
to be still excited, ai{t) is the probability amplitude for 
the transfer of the excitation to cavity mode oi, and so on 
as indicated by the defined photon creation and atomic 
raising operators. From this the density operator, at all 
times, is a sum of two parts: p{t) = po{t) + pi{t), with 


Po{t) = Pspon{t)\G){G\, 

(10) 

Pl{t) = |^(t))(^(t)|. 

(11) 


where the norm of |^^(t)) is 1 —Pspon(i), be., it is the null- 
record probability, the probability that no photon has 
been emitted. As C acts only in the one-quantum sub¬ 
space (the between-jump evolution), and V takes |'0(t)) 
to |G) (the quantum jump), we find from Eq. (|T|): 

pi(t)=Cpi, (12) 

po{t)='Dpi. (13) 

The equation of motion for the emission probability fol¬ 
lows from Eq. (US, 

=J22k. [|a.(t)|^ + m)?] +7a|6(04 

* (14) 

where Ki = KG-I-Kex , and Eq. (11211 yields the non-unitary 
Schrbdinger equation 

= (15) 

with non-Hermitian Hamiltonian 

H = Hi + H 2 + Hi2 — + J^Jb) 

- iKl^\a\ai + b\bi) - iK[l\ala2 + 4^2) 

+(^tcr2)- 


(16) 
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The nonunitary Schrodinger equation leads to a system of 
coupled differential equations satisfied by the probability 


amplitudes: 

6 = - (17a) 

di = - (ki + f Ai) ai - ig^^i - ihiPi, (17b) 

/3i = - (ki +?Ai)/3i 

- igi^i - ihlai - 2\J six'Kex /32, (17c) 

6 = - - 192^2, (17d) 

02 = — (k 2 + iA 2) a2 

- *526 - *^2/32 - 2 Y^( 17 e) 

^2 =- {k2 + *A2) P2 - 192^2 - ih2a2, ( 17 f) 


Emission spectra for the various outputs can be expressed 
in terms of the Laplace transform of these probability 
amplitudes. 


B. Emission spectra 

Spontaneous emission is a non-stationary process. Its 
spectrum can be computed from a double integration 
over a two-time correlation function, a formula derived 
from the response of a tunable oscillator driven by the 
emitted light [2^, e.g., a tunable filter cavity cascaded 
with the source of the spontaneous emission (Sec. 19.3.2 
of [2^). In the present case, the spontaneous emission is 
divided between several output channels; we treat them 
one by one. 


1. Side emission 

Emission out the sides of the microtoroids—to follow 
the nomenclature of cavity QED—consists first of atomic 
spontaneous emission to all modes other than ai and bi, 
i.e., atomic emission into free space. The atomic emission 
spectra are {i = 1,2) 

p(X> p(X> 

= dtj^ dt'e—(‘-‘')(a+(t)a-(t')), 

(18) 

where, from the quantum regression formula, with t' > t, 

= Tr [p{t)at] } . (19) 

More generally, the upper limit of the integral in Eq. m 
replaced by a finite time T and a spectral filter band¬ 
width, r > 0, may be introduced. Such time-dependent 


spectra for one-way cascaded inicrotoroids are calculated 
in [2^. 

Using the unraveling of Sec. IIII AI Eqs. (IT^ and (fTO 
lead to a spectrum expressed as the squared modulus of 
the Laplace transform of We may first write 

( 20 ) 

Noting then that {C+V) [|'!/i(t))(G|] = —iH\^|;{t)){G\, we 
find 


g(c+p)(t'-t) [|^(i))(G|] =e-*^(‘'-‘)|V;(t))(G| 

= m')){G\, (21) 


and hence, from Eq. (HU), 


(i')) = ?i'(0tr [cr, \ilj{t')){G\] 
= t'>t. 


( 22 ) 


With {af {t)<T~ (t')) = (a^ (0)* j result holds for 
t' < t as well. We arrive, thus, at the atomic spontaneous 
emission spectra {i = 1,2): 


rjiii) 

^ side,(7 


(w) 


7A 

27r 
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(23) 


where denotes the Laplace transform evaluated at 

s = 


nOO 

(24) 

Jo 

Spontaneous emission spectra for the four microtoroid 
scattering outputs might be defined and calculated in 
the same way. Most generally, however, the side emission 
might not be resolved and assigned to individual sources. 
The most natural generalization, built upon the method 
above, is therefore 


7side(^) — 


E 


|^7i'i|ii(w) -|- Aid'i{ui) -|- BiPi(ui) 


(25) 


with Ei, Ai, and Bi complex constants determined by 
the way the side-light is collected and hence interferes. 


2. Fiber emission 

The spectra of primary interest relate to emission from 
the microtoroid WGM’s into the fiber outputs. There are 
two fiber outputs, one, output a, traveling from left to 
right, and the other, output b, traveling from right to left. 
In parallel with Eq. m, the two fiber emission spectra 
are defined by 
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-I nOO /'OO 

Jfiber,a(w) = ^ y dt J dt''> Gout {t')) , (26a) 

-I POO POO 

'7flber,f,(w) = -^ J dt J dt'e~'‘‘^^*~*'\bl^^{t)bout{t')), (26b) 

where Oout and 6out are annihilation operators for output fields traveling in opposite directions [2^: 

aout(i) = -ain(t) + Ja{t), (27a) 

bontit) = -binit) + Jbit), (27b) 


where Oin and 6in are the corresponding inputs—in our case vacuum fields. By following the steps leading to Eqs. (1^51) 
and (1^ . the spectra are calculated as 


Efiber,a (^) 
d}iher,b(^) 



(28a) 

(28b) 


IV. STRONG-COUPLING REGIME 

We focus in this section and the next on two regimes 
of special interest in what, more generally, is a very large 
parameter space. We first consider strong coupling, i.e., 
conditions where the dipole coupling constants are much 
larger than the dissipative rates, where gi ^ {k^, 7 a}. To 
simplify matters, we assume the two microtoroid systems 
identical, with 

Ai = A 2 = A, Ki = K2 = K, (29) 

and also that intrinsic losses may be neglected compared 
against the coupling into the fiber, i.e., 

We thus write 

«: = K1,2 = (30) 

We also set 

hi = h 2 = 0, and 4>a = 4>b = 0. (31) 

With these simplihcations we can identify two distinct 
dynamical behaviors, each with its own distinct features 
written into the spectrum of spontaneous emission. The 
distinction arises from the relative phase of the dipole 
coupling constants, whether the coupling constants have 
the same phase or their phases differ by 7r/2. 


A. Coupling between delocalized atomic and 
photonic excitations: Qi — 02 — g 

In the first instance we assume equal dipole coupling 
constants, gi = g 2 = g, which (without loss of generality) 


we may assume real. We define 



= ;^ (6 ± 6) I 

(32a) 


T± = ^ [cii + [3i± (02 + P 2 )] , 

(32b) 


= 2 [“1 ~ PiT ( 0:2 — P 2 )], 

(32c) 

where for our assumed initial condition, Aj-(O) = l/-\/2 
and I±(0) = Z±i0) = 0. The equations of motion, 
Eqs. (|17al)-(ll7fp. decouple into a pair of independent sets 
of three: 


X+ = -^X+ - tV2gY+, 

(33a) 


= —2kYi. — iV2gX+ — kZ^, 

(33b) 


Z+ = kYj^, 

(33c) 

and 

X_ = -^X_-iV2gY., 

(34a) 


r_ = -iV2gX_ + kZ_, 

(34b) 


Z- = —2kZ- — kY-. 

(34c) 


The dynamics is dominated by the coupling at rate \/2g 
between X± and Y±. Notably this coupling is between 
delocalised atomic and photonic excitations— X± and Y± 
are sums or differences of excitation amplitudes applying 
to different microtoroids. 

From Eqs. (I28al) and (I28b|) . for the situation considered 
the spectra of spontaneous emission into the fiber may be 
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written in the form 


Tfiber,a(w) = - Y+{uj) + Z_{uj) 


Y+{UJ) - 


2k — iuj 


Y-i^) 


and 


(35a) 


Tfiber,6(w) = - V+(w) - Z-{uj) 


Y+{uj) + 


2k — iuj 


y_(cb) 


(35b) 


As shown by the example in the upper frame of Fig. [21 
these spectra are similar and dominated by vacuum Rabi 
resonances at frequencies close to a; — wa = ^\f2g. Note 
that the example takes smaller than Kex’ but not 

equal to zero, so spectra are computed from the full set 
of amplitude equations rather than Eqs. (I35al) and (I35bl) . 
Added to the vacuum Rabi doublet, the weak coupling of 
Y± to Z± give rise to small features around w — wa = 0. 
For clarity, these are shown on a logarithmic scale in the 
inset, where the broad feature (width ~ 2k) is associated 
with Z -, while the very narrow feature is associated with 
Z^ and a system eigenstate 


g{-a\ + b\+al- bl) - iK{a^ + ct^)] |G) 


^2^2 + 4g2 

(36) 

which is dark for emission into the fiber, i.e., Ja|'*A) = 
Jbl’fjj) = 0; with no emission into the fiber, spontaneous 
emission to free space modes occurs with a decay rate 
~ 7a(k/25)2 < 7A. 


take gi = g real without loss of generality. In this case 
it is instructive to consider the following combinations of 
excitation amplitudes: 


± + I3i), (37a) 

Y± =-^^2Ti^{a2 - P 2 ), (37b) 


and 


V = Pi), 

W = (q!2 + P 2 ) , 


(38a) 

(38b) 


where for our chosen initial condition, ^±(0) = l/-\/2 and 
F±(0) = 1/(0) = WpS) = 0. From Eqs. ([I^-dlZS), it is 
then possible to derive transparent equations of motion 
by invoking a rotating-wave approximation with respect 
to the large frequency g. We obtain, with the help of the 
said approximation, 






and 


i (^ + «)±*v^g 


X±+i^Y±, 

Y±+i^X±, 


(39a) 

(39b) 


1/ = -k1/ + kW, (40a) 

W = -kW - kV. (40b) 


B. Coupling between local atom-photon dressed 
states: g\ = ig 2 = g 

We turn now to dipole coupling constants that differ 
in phase by 7r/2, with gi = ig 2 = g, where again we may 

I 


In contrast to the previous case, the equations for X± 
and Y± now couple local normal-mode excitations at rate 
k/2; note the normal-mode frequencies iV^g as before. 
From Eqs. (I28al) and (I28bl) . the fiber output spectra are 
given by 


Tfiber.a(w) = — (w) -|-f 4 +(w) - X_ (w) -|-(w) 

47r L J L 

7fiber,b(w) = [-A+(w) - fE+(w)j - [x-(w) - iY-{u}) 

47r L J L 


V2 

V2 


W{uj) -k E(w) 


W{uj) t/(w) 


(41a) 

(41b) 


The coupling between normal modes manifests itself 
prominently in the spectrum of spontaneous emission 
into the fiber traveling from left to right, where each 
vacuum Rabi resonance takes the form of a doublet with 
splitting ^ K, as shown in lower frame of Fig.|21 Emission 


into the fiber traveling in the opposite direction shows 
no such splitting. The (small scale) structure around 
ui — uiA = 0 is determined primarily by the amplitude 
W + V. 
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FIG. 2: Spectrum of spontaneous emission into the fiber from 
the numerical solution of the full set of coupled amplitude 
equations; for (in units of 27r-MHz) 7 a = 5, pi = (p 2 ,*P 2 ) = 
50 (upper, lower), = 5, = 0.1, and all the other 

parameters zero. The blue solid (green dashed) curves plot 
the spectrum for the fiber output traveling from left to right 
(right to left). 


V. BAD-CAVITY REGIME 


The second regime of special interest is the bad-cavity 
regime, which corresponds to conditions where the field 
mode decay rates, ni and K2, are much larger than the 
dipole coupling strengths, gi and p2, and the rate of free- 
space atomic spontaneous emission ja- In this situation, 
one may adiabatically eliminate the field modes to arrive 
at a simpler model describing the atoms alone. We may 
do this starting from the coupled amplitude equations, 
Eqs. (I17a[) - (ll7fl) . or from the original master equation. 
We consider both approaches below. We first make the 
adiabatic elimination in the amplitude equations, from 
which we recover to two distinct types of effective atomic 
behavior: superradiant decay in the absence of dipole- 


dipole interaction, and dipole-dipole interaction without 
superradiant decay. We then show how this phenomenol¬ 
ogy emerges from the master equation after adiabatic 
elimination. The master equation is useful in its own 
right, as it facilitates our generalization from two to a 
string of microtoroids in Sec IVIl 


A. Adiabatic elimination of the cavity modes: 
amplitude equations 


To adiabatic eliminate the field mode amplitudes from 
Eqs. (ll7aD - (ll7fD . we set the time derivatives of ai, /3i, 
a 2 , and (32 to zero and solve the resulting set of algebraic 
equations to express the field mode amplitudes in terms 
of and ^2- Substituting the expressions into Eqs. (I17all 
and (|17dl) yields the pair of equations: 


• _ / 7A 2|gip ^ f 2 Kix niJ ^20*'^'’ 

^ V 2 Ki+iAiJ ^ {ki + iAi){K2 + iA2) 

(42a) 

and 

• _ / 7a 21^2 P ^ c 2 \J <7i <72 

^ V 2 K2 -I-IA2 / ^ (ki-I-1Ai)(k 2 + fA2) 

(42b) 


With the simplifying assumptions of Eqs. dsn])-®, and 
for the chosen initial condition, 5i(0) = 1 and ^2(0) = 0, 
the solutions are 


aw 

Ut) 

where 


with 


2 K + iA p 

2Kglg2 — e^-* 

{K + iA)'^ p 
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(43a) 

(43b) 


7A |gip + |g2p 

2 K + iA 2^’ 


(44) 


P 



\92\^f + 


f 2ti:g*g2 
y A? + zA 


2 


(45) 


For dipole coupling constants of the same magnitude, but 
possibly different phases, i.e., gi = e“®®g2 = g (real), the 
solutions reduce to the simpler form 


with 


aW = ^(e"+‘ + e"-‘), (46a) 

6 (t) = ^(e^+‘-e^-‘), (46b) 


lA _ 2g^ ^ 2Kg‘^e^ 

2 K + iA (K + iA)^^ 


( 47 ) 











































or, setting A = 0, 


0.012 


This expression for the dynamical eigenvalues reveals the 
two advertised regimes of control over the field-mediated 
atom-atom interaction, first superradiant decay in the 
absence of dipole-dipole interaction (0 = 0), and dipole- 
dipole interaction without superradiant decay (0 = 7r/2). 
More generally, from Eq. 621), the detuning can make the 
selection; for example, dipole-dipole interaction without 
superradiant decay is also selected with 0 = 0 and A = k. 
For simplicity, we consider only A = 0 in what follows. 


1. Superradiant decay: 0 = 0 


With both A and 0 set to zero, the dynamical eigen¬ 
values reduce to 


A+ — 


7A 
2 ’ 




(49) 


The term j k added to the eigenvalue A_ indicates 

two-atom superradiant decay [l^, (see also Sec. IIIA 
of [1^). A single atom coupled to one field mode in 
the bad cavity limit decays at the cavity-enhanced rate 
7 A/2-|-g^/K. Spontaneous decay through the two output 
channels of a single microtoroid (coupling to two counter- 
propagating modes) multiplies the enhanced decay rate 
by a factor of 2. The second factor of 2, resulting in a 
enhancement, comes from the collective character 
of the two-atom cavity-enhanced decay. The eigenvectors 
corresponding to the eigenvalues A± are 




^(cr+ ±cr+)|G). 


(50) 


These are the familiar subradiant and superradiant states 
of two-atom collective spontaneous emission. 

The spectra of spontaneous emission into the fiber are, 
from Eqs. (I28al) and (I28b|) with k^x = Kex = k and 

4^a — — 0 , 

7fiber.a(w) = -|di(w) -k a2{uj)\'^, (51a) 

TT 

2fiber,b(<j-') = — ■ (51b) 

TT 

Since, in a bad cavity with 0 = 0, 


ai + a2 = —(^1 + h) = — C 2 ), 

K 

both reduce to the Lorentzian 


^fiber,o(^) — -^fiber,6(^) 


ig^ 1 

TT K A?_ -f 


(52) 


Hence, spontaneous emission in both directions displays 
the superradiant linewidth for our assumed initial condi¬ 
tion. The upper frame of Fig. [3] compares the Lorentzian 


(a) 


0.008 

3 ; 

I 

bn 



( 0 - 0 )^ 

FIG. 3: Spectra of spontaneous emission into the fiber from 
the numerical solution of the full set of coupled amplitude 
equations; for (in units of 27r'MHz) 7 a = 5, gi = {g2, —192) = 
50 (upper, lower), = 500, = 0.5, and all other 

parameters zero. The blue solid (green dashed) curves plot 
the spectrum for the fiber output traveling from left to right 
(right to left). The red solid and dashed lines are plotted from 
Eq. ((52l) (upper) and Eqs. (I55al) and (I55bl) (lower). 


of Eq. ((5^ with numerical results computed from the 
full set of coupled amplitude equations. The agreement 
is very good, and improves with increasing k as expected. 
Note that the adopted parameters are similar to those of 
recent experiments with cesium atoms and microtoroidal 
resonators Q. 


2. Dipole-dipole interaction: 6 = 7r/2 


With A set to zero and 0 set to 7r/2, the dynamical 
eigenvalues become 


A± 


7A _ 2g^ ^ . 2g^ 
2 K K ' 


(53) 
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with corresponding eigenvectors given by Eq. dsni) once 
again. In this configuration, the cavity-enhanced decay 
is no longer superradiant, but replaced by a dipole-dipole 
coupling of magnitude 2\g \^/k between the atoms. With 
excitation amplitudes satisfying, 

ai+ a2= i — I3i+ 1^2 = —*-(Ci ~ *^ 2 ), 

K K 


the spectra of spontaneous emission into the fiber are 


^fiber,a(a^) 

Efiber,fe(a^) 


1 

9^ 

1 — i 

1 - 1 - i 

2 

(55a) 

47r 

n 

A_(_ “h iuj 

A_ -1- iw 

1 

1 

9 ^ 

1 - 1 - i 

1 — i 

2 

(55b) 

47r 

K 

A-(_ “h ioj 

A— -|- iuj 



These expressions yield different emission spectra in the 
two output directions, as illustrated in the lower frame 
of Fig. [3l For the chosen initial condition, only emission 
into the fiber traveling from left to right gives rise to a 
spectrum in the form of a doublet; with a peak-to-peak 


splitting of ^ the doublet is a clear indication 

of dipolc-dipole-induced frequency shifts. The spectrum 
of photon emission in the opposite direction exhibits a 
single broadened peak—the result of constructive rather 
than destructive interference of emission from the two 
collective modes. Note the comparison of the adiabatic 
elimination against the full set of amplitude equations, 
where the agreement is again very good. 


B. Adiabatic elimination of the cavity modes: 
master equation 


Finally, we use standard techniques (e.g.. Sec. 13.2.1 of 
[ 2 ^) to perform the adiabatic elimination in the master 
equation directly, where, working with the simplifying 
assumptions of Eqs. (ESI)-®, the master equation for 
the reduced density of the two atoms coupled through 
the adiabatically eliminated fields is found to be 


7A , 2k|5?| 


PA=-^[H,PA]+{f + ) D[a^W +[f + ) Dla^pA 


7A , 2K|g| 


- 2Re 


9291 


(k -|- iA)^ 


{2a^ PACft - <^2 PA - PAcrtcri + H.c.) , 


(56) 


with Hamiltonian 


H = -- 


2A 


: 2 -h A 2 


(I 51 I 0-1 CTi + I 52 I cr+cr 2 ) - 2Im 


929i 


{k + iA)2 


iafa2 + ^2 )i 


(57) 


The collective terms in Eqs. (|56)l and (|57l) [also (1621) and 
(ESI)] generalize the master equation for atoms interacting 
through a one-dimensional waveguide; in the latter case, 
the argument of the complex number gsPi ^/{k + iA)^ is 
replaced by ikA\zi — Z 2 \, with Pa the wavenumber at the 
atomic resonance frequency, and zi and Z 2 the atomic po¬ 
sitions [iglEljl . From the above general equation, master 
equations for collective decay without (coherent) dipole- 
dipole interaction and dipole-dipole interaction without 
collective decay are recovered as follows. 


The first term on the right-hand side describes two-atom 
superradiant decay. Note that for the dynamical eigen¬ 
vectors of Eq. EQl), 

J-|^+)=0, J-|^_) = y2|G). (60) 

Note also that the particular form of the collective jump 
operators in Eq. (1591) — with a minus sign—follows from 
choosing gi and g 2 to have the same rather than opposite 
phase (and the choice =</)& = 0). 


1. Superradiant decay: ff = 0 


2. Dipole-dipole interaction: 9 = 7r/2 


Setting both A and 9 to zero, with gi = g 2 = 9 (real), 
Eq. (I5S1) simplifies as 

PA = —D[J~]pa + ^ {D[a^]pA + D[a2]pA) , (58) 

Ki Zd 

where 

J±=cr±_^±^ (59) 


Alternatively, setting A to zero and 9 to 7r/2, Eq. (1551) 
simplifies as 


■V r 4 
PA = I -[crC 




(Jo -|- (Jo (7 


1 I Pa] 


+ ]pA + T’[ct2 ]pa)- (61) 
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This master equation describes cavity-enhanced, single¬ 
atom spontaneous emission in the presence of a cavity- 
mediated dipole-dipole interaction between the atoms. 


through a single optical fiber. Working again with the 
simplifying assumptions of Eqs. (l2^-(l3T]). Eqs. (I56l) and 
(EB generalize as 


VI. GENERALIZATION TO MANY SITES 
A. Master equation 

The master equation in the bad-cavity regime is readily 
generalized from two to N microtoroid systems coupled 

I 


N 


PA = 


= - i[H,PA] + 


k^l 

N N 

E E 

k^l 


7A , 


2 k‘^ + A^ 


D[<Jk ]PA 


9i9k 


l-k-l 


(k + zA)^ 


(2(7; pAcrt - CTfe PA - PACT+ a I + H.c.) 


with Hamiltonian 

H = 


- E E 

k k=l l=k+l 


2k 


9i9k 


J-k-l 


{k + iAy 




(62) 


(63) 


r 


where 


K — iA 
K + iA 


(64) 


As in the two-atom case, we readily identify interesting 
limiting cases for special choices of the dipole coupling 
strengths. We focus on the natural generalizations of the 
examples discussed in Secs. m and El 


each pair interaction shows a strong dependence on the 
relative positions of the atoms in the pair. The master 
equation may be written in the form 

N 

PA=- i[H, pa] + E '^^[^k\PA 

k^l 

+ ~^{^[Jod<^PA + D[Jeven\PA), (67) 


B. Special cases 

First, setting A to zero with gi = g 2 = ■ ■ ■ = g (real), 
Eq. (15^ simplifies to the master equation for iV-atom 
cavity-mediated superradiant decay: 

2 2 A! 

PA =—D[J-]pA + '^^D[ay]pA, (65) 

Hj z 

with 

J“ = (7f - (7^ -I-1- (-l)^+^(7^. (66) 

The second example from Secs. IIVI and El with a 7r/2 
phase difference between gi and 32 , does not generalize 
in quite such an obvious way. We generalize to alternat¬ 
ing real and imaginary coupling strengths: 51 = <73 = 
••• = —igi = —ige’’’ = 9 (real). This configuration 
no longer yields dipole-dipole interaction without super¬ 
radiant decay, but rather a combination of both, where 


with Hamiltonian 

E (68) 

and 

J^dd = +■■■ ^ (69a) 

Jeven = <^2 + ^4 + ^6 “I-■ (69b) 

Thus, each atom participates in superradiant decay with 
atoms located an even number of sites away, and, at the 
same time, dipole-dipole interaction with atoms located 
an odd number of sites away. 

C. N = 3 

As an illustration, we plot spectra for V = 3 in Figs. 4 
and 5. Results are computed from the full set of coupled 
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amplitude equations—i.e., the extension of Eqs. (I17al) - 
(fTTfl) to many sites—without the adiabatic elimination. 
The parameters, however, are chosen to approximate the 
adiabatic regime, and Eqs. (ESI) and EH) provides clear 
insight into the features observed. 

The spectra of Fig. 6 correspond to conditions where 
superradiant decay is predicted. They are dominated by 
a broad single peak of halfwidth ^ N2g‘^/K = 
with an additional small feature associated with single- 
atom spontaneous emission at rate 7 a. Notably, there 
is a dependence on the position of the initially excited 
atom, i.e., whether it is at the center of the group of 
three or at one end: the spectra depend on the direction 
of propagation when the excited atom is at one end, but 
are independent of direction—by symmetry—when the 
center atom is excited. Differences in the former case are 
seen from our numerical solution to become smaller the 
larger k becomes relative to the dipole coupling g. 

The spectra of Fig. [5] take gi = —ig 2 = 93 = 9, where, 
for three atoms in this configuration, the master equation 
in the adiabatic approximation, Eq. (EH), is 


PA=- i[H,pA] + ^{D[a^ ]+D[a2 ] + Dla^ ])pA, 

-h —D[crf -f ]pA + —D[(J 2 ]pa, (70) 

K K 

with dipole-dipole Hamiltonian 

2o2 

H = -(-cr ^^1 + 0-3 ) • (71) 


The atoms on the ends participate in cavity-mediated 
superradiant decay, while the center atom exhibits 
only cavity-enhanced single-atom spontaneous emission. 
While not participating in the superradiant decay, the 
center atom is nonetheless coupled to its neighbors 
through cavity-mediated dipole-dipole interactions. 

Diagonalizing dipole-dipole Hamiltonian dZD yields 
the pair of eigenkets 

|e±) = T V2cr+ - cr+) |G), (72a) 


with eigenvalues e± = ±2^g^ j k, and the eigenket 

|eo) = ^ (a+ + 4) |G), (72b) 

with eigenvalue eo = 0. Features appear at correspond¬ 
ing frequencies in the spectra of Fig. [S] Note the action 
of the jump operator +a^ on the three eigen¬ 
kets: JJddle-t) = = 0 and J”dd|eo) = V^jG). 

Thus, with initial state cr+|G) = (|e+)-|-|e_)-|--\/2|eo))/2, 
one observes contributions to the emission from all three 
eigenstates, albeit in different combinations depending 
on the emission direction. On the other hand, for initial 
state cr^lG) = (l/-\/2) (|e+) — |e_)), superradiant decay 
is absent from the dynamics and the spectrum takes a 
straightforward double-peaked form as a result of the 
dipole-dipole interaction between neighboring atoms. 



(jO-CjO, 

A 



A 

FIG. 4: Spectra of spontaneous emission into the fiber in 
the bad cavity regime for three atoms. Curves plotted from 
the numerical solution of the full set of coupled amplitude 
equations; for (in units of 27r-MHz) 7 a = 5, gi = g 2 ~ 93 = 
50, = 500, = 0.5, and all other parameters 

zero. The blue solid (green dashed) curves plot the spectrum 
for the fiber output traveling from left to right (right to left). 
The initial state is o-^\G) (upper) and crJ|G) (lower). 


VII. CONCLUSION 


We have taken a first step towards characterising and 
manipulating the emission properties of two-level atoms 
through their coupling, via evanescent fields, to a net¬ 
work of fiber-linked microtoroidal resonators. Using the 
theory of two-way cascaded systems, we computed the 
spectrum of spontaneous emission into the fiber outputs 
for a pair of atom-microtoroid systems and illucidated 
its dependence on the relative phase of the atom-cavity 
coupling: under strong-coupling conditions, the decay is 
dominated by an interplay of either delocalized or lo¬ 
cal excitations, depending on relative phase, while under 
bad-cavity conditions (microtoroid and fiber fields adia- 
batically eliminated), the atoms exhibit either collective 
spontaneous emission with no dipole-dipole interaction, 
or a (coherent) dipole-dipole interaction and single-atom, 
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(O-CO, 

A 



( 0 - 0 ), 

A 


FIG. 5: As in Fig. |4]but with gi = —ig 2 = (?3 = 50. 


cavity-enhanced spontaneous emission. In the latter case, 
we recover behavior expected for atoms coupled through 
a one-dimensional waveguide from our two-way cascaded 
system model. 


We extended our results under bad-cavity conditions 
to many two-way cascaded atom-microtoroid systems. 
While the conditions leading to collective spontaneous 
emission with no dipole-dipole interaction carry over, a 
more complex combination of collective decay—involving 
atoms separated by an even number of sites—and dipole- 
dipole interaction—atoms separated by an odd number 
of sites—is seen in the second case. 


An externally driven network raises interesting issues 
untouched in this work; it may attain a non-trivial steady 
state, e.g., one showing significant entanglement between 
the atoms Preliminary results on the driven 

version of our two-way cascaded system are presented 
in [l^ . though this is a direction wide open for further 
study. 
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